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Abstract 



We consider the computation of the entanglement entropy in curved back- 
grounds with event horizons. We use a Hamiltonian approach to the prob- 
lem and perform numerical computations on a spherical lattice of spacing 
a. We study the cosmological case and make explicit computations for the 
Friedmann-Robertson- Walker universe. Our results for a massless, minimally 
coupled scalar field can be summarized by S en t = 0.30r^/a 2 , which resem- 
bles the flat space formula, although here the horizon radius, th-, is time- 
dependent. 
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I. INTRODUCTION 



The formal proof of the validity of the first law of black hole thermodynamics [|1],0] for 
arbitrary perturbations of a stationary black hole in any gravitational theory derivable from 
a diffeomorphism invariant Lagrangean gave an impulse and renewed interest on the com- 
putation of the black hole entropy. Recently, several authors have obtained expressions for 
the thermodynamic or coarse-grained black hole entropy in the form of an integral of the 
entropy density on any cross-section of the Killing horizon. More explicit expressions have 
been given for different higher derivative gravitational (interacting with matter) Lagrangeans 
0-0. The main techniques used to derive the black hole entropy have been the Noether 
charge construction ||, the Euclidean signature method [0, and the field redefinition ap- 
proach 0. 

All the above derivations make use of the first law of black hole thermodynamics where 
all its terms are known (mass, temperature, angular momentum, angular velocity of the 
horizon, electromagnetic charge and potential), but an integral on the horizon. One pre- 
cisely identifies this unknown quantity on the horizon with the black hole entropy by analogy 
to any other ordinary thermodynamic system. This can be called the coarse-grained en- 
tropy because it is derived in the thermodynamic limit. It is, however, desirable to have 
an independent, statistical mechanical derivation of the entropy, which may be called the 
fine-grained entropy. Relatively early attempts PJToH to compute this quantum-originated 
entropy considered the statistical entropy, 5* ~ k B J2nPn^Pn ~ k B \iaN, with N being the 
number of modes or configurations of field very close to the black hole horizon, 't Hooft H 
has computed the quantum corrections to the black hole entropy by counting the number 
of classical eigenmodes of a scalar field outside the event horizon. He immediately noted 
the appearance of an ultraviolet divergence in the entropy as one approaches the horizon in 
counting the modes, and discussed the need of imposing a short distance cut-off. 

Bombelli et al. |TT| studied instead the entanglement entropy S = — Tr [p In p] by con- 



sidering in flat space a scalar quantum field in its ground state and computing its reduced 
density matrix, p = |$ >< $|, where one has traced over the field degrees of freedom 
inside a region of space that is intended to represent the interior of a black hole. Thus, 
the entanglement entropy is associated to correlations across the boundary separating two 
regions. More recently, several fl2HrTf computations of the entanglement entropy agree in 
obtaining the proportionality to the area A of the surface dividing the known from the un- 
known regions and the 1/h 2 ultraviolet divergence, where h is the minimal distance to the 
horizon we reach in counting modes of the external field, typically given in Minkowskian or 
Schwarzschild radial coordinates. Their results can be summarized by 

where g accounts for the number of fields and its helicities and 7 = l/(3607r) (see Refs. 
0,0)- 



In Refs. p^JT7UT8[1 the entropy is computed by identifying the dynamical degrees of 



freedom with the states of quantum fields propagating in the black hole's interior. In Refs. 
the entropy of fields outside the black hole is studied instead. The two kinds of 



computation are conciliated in Ref. JT^| where it is shown that the density matrices of the 
two situations are identical. 
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The question of the ultraviolet divergence of the entropy is discussed in Ref. |TJ| where 
the Pauli-Villars regularization is used to avoid the introduction of a Planck scale cut-off in 
the computation of the entropy. In the paper ]2(J the entropy difference between two states 
is shown to be better behaved. Susskind |2l| speculates about the possibility that the black 
hole entropy be finite in string theory. 

In the present paper, we handle this problem, following Ref. |13| , by making computations 
on a lattice of spacing a, what renders the entropy finite. In Srednicki's work the reduced 
ground state density matrix of a free, massless, scalar quantum field is obtained by tracing 
over the degrees of freedom located inside of an imaginary sphere of radius R and by making 
use of the analogy with the result 1TJJ for a coupled system of harmonic oscillators. The 



final result is 



S 



0.30 {4nR 2 



47T 



(2) 



which, as the author remarks, bears an striking similarity with the formula for the coarse- 
grained entropy of black holes 



S 



A 



H 



A/2 ' 



(3) 



where Ah is the area of the black hole's event horizon. The similarity between these two 
equations is in fact surprising, since Srednicki's computation has been carried out in flat 
spacetime while we know that the geometry generated by a black hole is curved. The 
argument here is that the modes which contribute mostly to the entropy are those of short 
wavelengths that do not "feel" whether the spacetime is curved ||11|| . However, the event 
horizon is a null surface which greatly differs from the spacelike surface R = constant. In 
addition, we expect the entanglement entropy of a quantum field computed in the curved 
metric of a black hole to give just a correction (due to quantum fluctuations) to the coarse- 
grained entropy produced by the classical curved black hole background, but not the whole 
value (given by Eq. (||).) 

Comparatively much less attention has been paid to the cosmological case. Here one 
can also define a horizon temperature and the laws of thermodynamics [^2,23], although 
their interpretation is not as transparent as in the black hole case. In our paper we will 
concentrate on the entanglement correction to the classical entropy of the universe. In Sec. 
II we generalize the approach of Ref. [13| to the case when the spacetime is curved and 
possesses a horizon acting as a barrier for our knowledge of the quantum state beyond it. In 
Sec. Ill we review some properties of the Friedmann-Robertson- Walker geometry while in 
Sec. IV we give our main results. The paper finishes with a discussion of the cosmological 
results and it possible extension to the black hole case. 



II. COMPUTATION OF ENTANGLEMENT ENTROPY IN CURVED 

SPACETIME 

In this paper, we want to extend the existing schemes for calculating the entanglement 
entropy to a curved spacetime with horizon. To compute the entanglement entropy, one 
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considers a quantum field in its ground state. The spacetime is divided into two parts which 
are separated by a boundary (in spacetimes with horizon, this subdivision appears quite 
naturally). A given observer has access to the information of only one of these two regions. 
In his/her subsystem, the field cannot be described by a pure state wave function. Because 
of the correlations that exists between the two regions, even in the ground state of the field, 
it must be described by a density matrix with an associated nonzero entropy. This is the 
origin of the entanglement entropy. 

In Refs [0,0 a scheme was given for calculating the entanglement entropy in Minkowski 
space. They considered the quantum field as a set of coupled harmonic oscillators with 
Hamiltonian given by 



(4) 



In the usual flat spacetime there is not horizon and the boundary can be chosen arbitrarily. 
From @, it is then possible to compute the entanglement entropy numerically |13| . 
We consider a curved spacetime with an isotropic metric of the following form 



ds 2 = C 2 {r, t)dt 2 - B 2 {r, t){dx 2 + dy 2 + dz 2 ). 
For a massive, non minimally coupled scalar field 0, the Lagrangean is given by 



C 



-fj 



f(^)(^)-(m 2 + (i?)^ 



which becomes for the metric 

1 



-CB 



C- 2 (d t <P) 2 - £- 2 (V0) 2 - (m 2 + ^R)0 2 



The canonical momentum conjugate to <fi is 

dC 



n 



gg^d, 



C^B'^r 



<9(<9 O 0) 

We obtain the Hamiltonian of the scalar field from the definition 

H = J d 3 x [Ud t (j) - C] , 

which becomes in our case 



H 



d 3 x 



CB- 3 U 2 + CB(V<p) 2 + CB\m 2 + £R)4>* 



(5) 



(6) 



(7) 



(8) 



(9) 



(10) 



The problem can be reduced to an effectively one-dimensional one if we exploit its spheri- 
cal symmetry. The metric (|5]) has an isotropic form then decomposition in terms of spherical 
harmonics is possible. Following Srednicki in using the real spherical harmonics Zi m (8, ip) 
which are orthonormal and complete. We decompose the field and the canonical momentum 
into partial wave components 
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x 



I III 



n(x) = -Y,z lm {e,<p)u lm {r). 



1 1 a 



When we substitute ([Tl|) and (|12|) into (0), we obtain after some algebra 

/ o , \ 2 



2 ^ 



dr 



czr 3 nL + 



-'Zm 



6>T 



(11) 

(12) 



(13) 



We next perform a substitution in the canonical variables in such a way that the commutation 
relations are preserved: 



n; 



i in 



FIT 



Imi 



lm 



Hmi 



(14) 



with F(r,t) = C l l 2 B 3 / 2 . In terms of the variables (j)' lm and U.' lrn , our Hamiltonian takes the 
form of Eq. (W) 



h = -Y. 

2 ^ 

lm 



dr 



U'l+CB[—C^B 



d 



lm 



dr 



'lm 



+ 



C_ (d_ 
B~z [dr 



c 2 i(i + 1; 

.B 2 r 2 



+ c 2 K + wL 



(15) 



We can write this as 

2 ^ 

im 



2 (71+74) ,/ 2 , 72 ^ 9 A/ 
ll lm ' o rim ~r 



+73 



-7m 



<9r 



+ 75— ~ 2 (t>lm +760« m 



(16) 



where 



7i(r) = a 2 C5 (^B 



72 (r) = 2aC^-i£- (C^F 



id_ 
dr 



2 d-2 



7s(r) = 75 (r) = C 2 £ 



74(H 



«L C B-*£-(CB), 76 (r) = C 2 a 2 (m 2 + £J2). (17) 



In Eqs (0) and (0)) we have introduced an arbitrary length scale a to make the ji dimen- 
sionless. In the form (|T6|), the Hamiltonian is in a suitable form for discretization. 

In order to compute the entanglement entropy numerically, we put the system on a 
spherical lattice with iV points in radial direction. The lattice spacing is denoted a so that 
r = ja, where j is an integer. Discretization of ([TB]) leads to 



h = — y 

2a V 

jlm 



n Lj + ( 7i + 74 + 7s 



1(1 + 1) 

o'2 



+ 76 



Imj 
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+l24>lmj((f ) lmj+l - 4>lmj) + lz{.<Plmj+l ~ 4>lmjY 



(18) 



To discretize the parameters 7$, it is necessary to specify the form of the spacetime metric. 
This will be considered in the next sections. The Hamiltonian (|TJ) is of the general form 
(|) with 



K 



7i (j) ~ 72O') + 73 (j + 1) + 7a(j) + 74(j) + 7s (j) 



1(1 + 1) 



+ 7e(j), 



(19) 



Knowing the explicit form of the 7$, it is thus immediately possible to calculate numeri- 
cally the entanglement entropy using Srednicki's algorithm |13|]: The lattice is divided into 
two parts. We want to trace over the first n lattice sites (the region inside the horizon). To 

this end, the matrix K is written via singular value decomposition as K = U T KpU where 

1 

Kd is diagonal and U orthogonal. One defines the square root of K as u = U T K^U and 
writes it in the form 




(20) 



\B T A- X B 



where the submatrix A is nxn and C is (N—n) x (N—n). Next we introduce B 

and 7 = C — B. We perform a singular value decomposition of 7: 7 = V t ^dV- We finally 

_i _ _i 

define the matrix B' = , ~f D 2 VBV Tr Y D ' 2 which has N — n eigenvalues B[. The entanglement 
entropy associated with the trace over the n sites inside the horizon is 



£ = £(2/ + 1)3, 
1=0 



where 



N—n 

i=i 



log(l-6 



i-e, 



log & 



^.2 



(21) 



(22) 



and £i = B[/{1 + y 1 — B' A. A derivation of this algorithm as well as its physical discussion 



can be found in 13 



III. GEOMETRY OF ROBERTSON- WALKER UNIVERSES 



It is our aim in this paper to calculate the entanglement entropy of a quantum field 
in a curved spacetime with horizon. Specifically, we are interested in the behavior of the 
entropy for different classes of cosmological models. We consider therefore the Friedmann- 
Robert son- Walker (FRW) line element 



ds 2 = dt 2 - l3 2 C 2 {t) 



dx 2 



sin 2 (v / ^x) 



(d9 2 + sin 2 9dip 2 ) 



(23) 
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and treat the three different cases k = 1, —1, corresponding to a closed, open and spatially 
flat universe. The constant (3 has the dimension of a length and determines the spacetime 
curvature. 



In the form (|2~3"D , the metric is not very useful for calculating the entanglement entropy 
because it is not in the isotropic form (El). It can obtained by means of the coordinate 



transformation 24 



drj = C~ l (t)dt, 



x = —^L tan (Vk—^\ sin 9 cos (p. 
Vk V 2J 



(24) 



y 



Vk 



tan ( Vk-) sin^sinyj, 



Vk V 2,' 



Note that it does not mix space and time variables. For the sake of simplicity, we restrict 
ourselves in the following to times r\ > 0. 
With (|24T), the line element becomes 



ds 2 = C 2 {r)) 



drj'' 



r(dx 2 + dy 2 + dz 2 



(25) 



It is of the form (pj) and it will be the starting point for our subsequent calculation of the 
entanglement entropy. 

An important concept for our discussion is the notion of the cosmological event horizon. 
It represents the maximal coordinate distance a photon emitted at time t can travel. For a 



metric of the form (23), it is given by [25 



Vk 



2(3 



where 



At? = Tj{t f ) - v(t) = J 



if df 



(26) 



(27) 



Here, tf equals infinity for k — 0, —1, while for k = +1 it corresponds to the time when the 
universe collapses to a final singularity. 

We will be interested in the area spanned by the cosmological event horizon. Generally, 
in a curved spacetime with diagonal metric, the proper area of a spherical surface r = ro = 
constant is defined by [^6 



An 



For the metric (E5f), this becomes 



r=r 



V9ee9 w d6d(f. 



(28) 



A = 4nC 2 (t f n2{ ^ Xo) =4nC 2 ( V ) 
k 



2 \2- 



(29) 
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As a concrete example, one can consider the evolution of a universe with negligible 
matter contribution but endowed with a cosmological constant A. In that case, we have for 
k = ±1,0 the following evolution 



C(t) = \ (eW + ke-^) 



(30) 



where /3 — y 3/A. In terms of the conformal time 

Ar? = —^= arctan \pke~ 1 ^ 
Vk 1 

we have 



tan 



The "radius" of the cosmological horizon can then be explicitly computed from (E§): 



k 



C - VC 2 - k 



The area of the horizon can be calculated by inserting (|33"D into 

A H = 4tt(3 2 . 



We obtain 



(31) 



(32) 



(33) 



(34) 



It is noteworthy that it does not depend on C(r]), i. e. it is time independent. 

In the k = case, the evolution law is C(t) = exp(t//?)/2, corresponding to C(rj) = (3/ Arj. 
The horizon radius is 



r H = Ar]= —. 



(35) 



and the area is again given by (|34T). 



IV. ENTANGLEMENT ENTROPY IN FRIEDMANN-ROBERTSON- WALKER 

UNIVERSES 

With the formalism described in Sec. II and the geometry reviewed in Sec. Ill, it is 
now possible to calculate numerically the entanglement entropy of the quantum field which 
is associated with the FRW cosmological horizon. Let us first note that the metric ( |25| ) has 
the form @ with 

C = C( V ), B = -^L. (36) 

Evaluating Eqs. (|ITD, we find that the coefficients ji are 

9k 2 j 2 3kj{kj 2 + 4b 2 ) 



S 




■ , , , , kkj 2 + Ab 2 _, 2 2 

73 = 7s = ( 1 + k— } , 74 = - A , 7 6 = Ca 



m 2 -^(C + kC) 



■ (37) 



Here, the dimensionless quantity b is defined by b — (3 /a. Note that in the case m = = £ 
the 7i are independent of rj, as will also be the Hamiltonian. We shall now consider this 
simpler case. 

In Minkowski space, the radius up to which the field has to be traced over has been 
chosen arbitrarily jT3]. The entropy was then calculated for different values of this radius. 
In a spacetime with horizon, however, this procedure is not appropriate. The horizon sets 
a natural radius for the region to be traced over since it represents a boundary for the 
information an observer inside (in the cosmological case, outside for a black hole) can have 
about the outside region (inside for a black hole). We therefore proceed in the following 
way: At a fixed value of C(rj), we consider a set of universes with different values of (5. 
According to fl33|) , they possess different horizon radii r#. We calculate the entanglement 
entropy numerically for each of these universes by tracing up to the horizon radius fl3"B| ) (this 
is equivalent to tracing for radius bigger than the horizon due to the validity of the identity 
of the density matrices in both cases \TE\). For a special value of the parameter C, the result 



is plotted in Fig. | as a function of the proper area of the horizon which depends on f3 via 
([34]) . We find straight lines in all of the cases k = ±1, 0. This shows that the entanglement 
entropy is proportional to the horizon area for the three spacetime classes: 

S = M0£p (38) 

The slope of these curves is given by the a priori unknown quantity fk(C)/ (Atto 2 ) . It depends 
on k as well as on the expansion parameter C(rj). We can try to find the functional form 
of fk(C) by repeating the above procedure for different values of C and plotting f k versus 
C. Doing so, one finds the results shown in Fig. ||] where the dots represent our numerical 
computation of f(C) for k = ±1,0. 

It would be desirable to know the analytical form of the functions fk{C) and therefore 
of the entanglement entropy. Within our numerical approach, we can only try to find a 
function that matches our data points sufficiently well. In Fig. |, the solid lines are plots of 
the functions 

h(C) = 0.30^f = 0.30 • ^(C - V^k) 2 . (39) 

We see that the numerical values are perfectly fitted within the numerical precision. Note 
that there are no free fitting parameters apart from the numerical factor 0.30. For k = 
we recover the flat space Hamiltonian and f{C) = 0.30C -2 . The first equation in (p9[) 
holds for general FRW spaces, while the second equation, with the particular dependence 
on C(r]), corresponds to the simplified case where we can neglect the matter contribution to 
the evolution of the universe (de Sitter-like geometry). 
Taking into account (|3"3"| ) and (0), this can be written 

2 

S ent = 0.30^f . (40) 
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This simple formula is our final result for the open, asymptotically open and closed FRW 
universe. Formally, it resembles very much the Minkowski space formula @, although one 
should keep in mind the conceptual differences due to the curvature of the spacetime and the 
existence of an event horizon, which are also reflected in its specific calculation, as explained 
above. 



V. DISCUSSION 

As we have seen, in the conformal time gauge for the metric, Eq (p5|) , the Hamiltonian of 
a massless, minimally coupled scalar field is time independent. However, since the boundary 
(th) up to which we trace over depends on time, the entanglement entropy varies as the 
universe evolves. In particular, it will increase as the universe grows and decrease as it 
shrinks (for k — 1) in the same way the boundary of the observer's knowledge (given by the 
cosmological event horizon) does. The dependence on the particular way we take the slices 
of time has its origin in the Hamiltonian approach we have used. 

It is also interesting to recall here that in Ref || it was found, for black holes in quadratic 
theories of gravitation, also a proportionality between the entropy and the square of the event 
horizon radius rather than to the horizon area. Thus, in some sense, Eq (EOT) seems to be 



a more fundamental relation than Eq ([I]). Note however, that the above proportionality 
for black holes in quadratic theories holds in the conformal approximation to the solution 
made in Ref ||, but it may well not survive in an exact solution (see Refs @,[J). The other 
possibility to further investigate this point is to directly apply the procedure developed in 
this paper to the black hole case. This is presently under study by the authors. 
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FIGURES 



FIG. 1. This figure shows the proportionality of our computed entanglement entropy to the 
cosmological horizon area in a given slice of time, i.e. for fixed conformal factor C(rj). Here 
k = 1, —1,0 correspond respectively to a FRW closed, open and spatially flat universe. 



FIG. 2. Here we plot the function /fc(C) appearing in Eq (38) and defined by the proportionality 
entropy-horizon area. The dots correspond to our numerical computation while the continuous 
line to the fit by 0.30r^//? 2 . The k parameter reflects the effect of different spatial curvatures. The 
dependence on C(rj) reflects time dependent character of the de Sitter metric. 
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